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Abstract. We consider the energy spectrum of the planar quantum well which consisted 
of two ribbons of narrow-gap semiconductors and a graphene ribbon between ones. It is 
shown that the gapless mode appears only in case of inverted narrow-gap semiconductors. 
Spin splitting of the energy spectrum for a nonsymmetric quantum well is calculated 
taking into account a specificity of graphene. We investigate interface states and optical 
transitions. It is shown that the optical transitions are possible only with a conservation 
of a parity. 
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After the first experimental investigation of graphene, a monatomic layer of carbon atoms 
forming a regular hexagonal lattice [1, 2], an intensive activity has arisen in different directions. 
In particular, a multitude of works on electronic properties of narrow graphene ribbons with 
nanometer sizes (nanoribbons) of late three years were made. Among early ones, there were 
works on study of electron states of the graphene ribbons using the Dirac equation with the 
appropriate boundary conditions [3, 4]. The electronic properties of the graphene nanorib- 
bon (GNR) depend strongly on size and geometry, as well as on whether hydrogen atoms 
were deposited on free atomic orbitals of carbon atoms at edges of the GNR [5]. In their 
transport properties, the GNRs are very similar to carbon nanotubes, since a free motion is 
one-dimensional (ID) in both cases [6, 7]. 

The boundary conditions for nanotubes, in which charge carriers are described by the Dirac 
equation, were discussed in detail in Ref . [8] . Two types of the boundary conditions correspond 
to two types of edges of the graphene ribbon, namely zigzag and armchair [3]. A formation 
of the zero-energy surface states is characteristic for the GNRs with the former type of edges. 
The energy spectrum is gapless (metallic) for the latter type of edges, when the ribbon width 
d = (3M + l)a, where M is an integer and a = \/3ao is the lattice constant (ao is the interatomic 
distance), and insulating otherwise [4]. 

Recently, a transistor was made on the basis the GNR with the 2-nm width and 236-nm 
length (nanoribbons with the 10-60- nm width were also investigated) [9]. The GNR was taken 
so narrow to guarantee a wide semiconducting energy gap for room temperature operation of 
the transistor. However, it loses slightly a compactness to the transistor based on a graphene 
quantum dot with the 30-nm diameter [10]. 
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A planar contact of the GNR with zigzag edges with metal ribbons has been investigat- 
ed in Ref. [11] within the tight-binding approximation. In our paper, we consider the size 
quantization of massless charge carriers in the GNR having armchair edges with the width 
d 7^ (3M + l)a between narrow-gap semiconductor ribbons. Changing narrow-gap semicon- 
ductors, we can change potential barriers heights which are the band gaps of the narrow-gap 
semiconductors, what opens additional possibilities for the energy band-gap engineering [12]. 

The potential barrier height for the single GNR has been calculated in Ref. [13] within 
the local spin density approximation. It equals 2.33 and 3.33 eV for the armchair and zigzag 
edges, respectively. These values is one order higher than a typical band gap of the narrow-gap 
semiconductors. 

We carried out calculations within the Dirac model [14]. Let us direct X axis parallel to 
interfaces of the narrow-gap semiconductors and graphene, Y axis transversely to a heterostruc- 
ture plane, and Z axis transversely to X and Y axes (Fig. 1). In a two-band approximation, 
a heterostructure composed of narrow-gap semiconductors 2 and graphene is described by the 
Dirac equation. Half the band gap A, the work function V, and the matrix element for the rate 
of interband transitions u are constant and change only on the heterostructure boundaries. For 
chosen orientation of axes in given case, the Dirac equation has following form 

h d v = {uaV% + u a °^% + 7% + v<} * = ev, (1) 

where 7 °H- " laod-rH ° ') are t ee Dirae prices, J , t ee 2 x 2 unit mate , 




(T 

a are the Pauli matrices, p = — «V is the momentum operator, and h = 1. 
It is simply verified that the operator 3 



P = i 7 ° 7 V (2) 

commutes with Hamiltonian Hp. 

To clear up the sense of the operator P we present it in the form 

P = *7 A„ (3) 

where «7° is the inversion operator, = e~ 4 ^ Sy = — %L y is the operator of the rotation about 

the Y axis by angle 7r (see, e.g., [17]), and H y — I y ) . From this, one can see that the 

operator P plays the role of the pseudoparity operator introduced for the layered narrow-gap 
semiconductor heterostructures [18]. 



2 We can use a finite-gap graphene modification as a narrow-gap semiconductor [15]. 

3 One should distinguish the operator P from the helicity operator h, introduced for unbounded graphene 
[16], which is the operator of a projection of pseudospin onto a free motion momentum. The free motion is ID 
in the planar heterostructure and the ID operator h (written in the 4x4 matrix form) does not commute with 
Ho- Therefore, an eigenvalue of h does not conserve. The conserved quantum number is the eigenvalue of the 
operator P. 
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Fig. 1. Considered hetero structure: the GNR between two the narrow-gap semiconductor 
ribbons (here Y axis is directed from us; an atomic structure of narrow-gap semiconductors is 

not shown). 

The eigenvalues A = ±1 of the operator P determine sign of a spin projection onto Y axis. 
It can be verified by the calculation of the mean value of spin as it was made in Ref. [18]. 
Hereinafter, we will refer to A as the parity. As it was emphasized in our previous paper [19], 
the charge carriers in graphene are the chirality massless Dirac fermions (spin of particle is 
aligned with its momentum or is opposed to it). As it will be shown below, there is no the 
gapless mode in the quantum well considered. The charge carriers in the GNR of a finite width 
bounded by layers of usual (noninverted) narrow-gap semiconductors, have nonzero effective 
mass and its energy spectrum is similar to the semiconductor spectrum with the finite band 
gap. 

We find a solution to equation (1) in a class of functions being the eigenf unctions of the 
pseudoparity operator (2). Components of the eigen bispinor interconnect by relations 

^A2=«A^A1, 

(4) 

?/>A4 = -«A^A3- 

Now the Dirac equation (1) is written in the form 

{m;7°7 3 R - i\l 3 Uik x + 7% + V*} * A = E x ^ x . (5) 
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Two independent components of the eigen bispinor satisfy the system of equations 



~ iUi Itz + i ^ Ui ^ x ) = ( Ex + Ai ~ V ^ ^ A3 ' 

-ilk j ~ i\Uik x ) t ' A :; = ( t\\ - A, - r, ) ( . 



(6) 



We find the solution in three ranges in the form 

1) z < 0: Vai = A ie fclZ , 7/> A3 = A 3 e^, 

2) 0<z<d: Vai = #ie ifcz + £ x e- ite , ^as = 5 3 e ifcz + i^e"^ 2 , 

3) z > d: Vai = C ie -*»*, V>A3 = C 3 e~ k ^, 

respectively, for the narrow-gap semiconductor with Ai, V±, ui, graphene with A2 = 0, V2 = 
0, 1*2, and the narrow-gap semiconductor with A3, V3, M3. 

Using the boundary conditions of continuity of \/u^/ [20], one can exclude constants 
Ai, A 3 , Ci, C 3 , -B 3 , _B 3 and obtain the system of equations on coefficients B 1 and Z?i 1 



112k + i\u 2 k x iu\k\ — %\u\k 

E~x + 
U2k + i\v,2k x 



E x + A 1 - V 1 
iu 3 k 3 + i\u 3 k x 
E\ + A3 — V 3 



Bi4 



— u 2 k + i\u 2 k x iu\k\ — %\u\k, 



E x ' Ex + A.-V, 

U2k + i\u2k x iu 3 k 3 + i\u 3 k_ 



£i = 0, 



V Ex Ex + A 3 -V 3 J 

(7) 

The dependence of energy on the allowed wave number k and the wave number of the free 
motion k x is given by relation 

E x = ±u 2y /k* + kl (8) 

The gapless mode corresponds to k — 0. Then the energy spectrum is linear one E\ = ±U2\k x \. 
However, from (7), one can see that then Bi + B\ = and B 3 + B 3 = 0, what means that 
^x = 0, therefore this solution is not physical one. 

The gapless mode in the quantum well based on graphene is possible if graphene is bounded 
by the inverted narrow-gap semiconductors [21] in which AxA 3 < 0, u% = U2 = u 3 = u and 
Vi — V2 — V 3 — 0. The wave function is constant in the range 0<z<datk = 0, namely, 
ip xi = Di, tpx3 = D 3 . After a substitution of this solution in (6), we obtain D 3 = ±iDi and 
Ex = ±\u\k x \. Joining solutions on boundaries z = and z = d, we find the wave function 



1) *f ) 

2) *f ) 



3) ¥ x ±] = C¥ x ±] e^( 



z-d) 



where 



(±) 



/ 1 \ 

iX 
±i 
\±XJ 



4 The system of equations in B 3 and B 3 is analogous to (7). 
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A constant C is determined from the normalizing condition 

dz¥ x V x = 1, (9) 



u ( 1 1 



1-1/2 



d 



.2 VIAil |A 3 | 

The case D 3 = iD\ corresponds to Ai < 0, A 3 > 0, and E\ = Xuk x . For D 3 = —iDi, we obtain 
Ai > 0, A 3 < 0, and E x = -Xuk x . 

We consider below the case of the noninverted narrow-gap semiconductors with A 1.3 > 
when there no the gapless mode. It is simply obtained the equation in k corresponded to the 
equation investigated in Refs. [22, 23] 

tan k d = MC +5 } ^ (10) 

1 + B(C -C)-CC 



where 



A - — B - ^ U2 ^ x C - Ulkl - ^ Ulkx C " 3 ^ 3 + ^ U3 ^ 



E x ' E\ Ex + Ai-VS E x + A 3 -V 3 ' 



u\k\ = A? - (E x - VJ 2 + ulk 2 x , u\k 2 = (E x - V 2 ) 2 -A 2 2 - ujk 2 x , u 2 3 k 2 3 = A 2 3 - (E x - V 3 ) 2 + u\k 2 x . 

In general case, equation (10) is sufficient awkward therefore we restrict oneself to a consid- 
eration of the particular cases for the defined types of the quantum wells. 
The symmetric quantum well. 

We refer to such quantum well in which V\ = V 3 = V, Ai = A 3 = A, A 2 = 0, V% — 0, U\ = 
u 3 = u, 112 — u as the symmetric quantum well. Equation (10) is reduced to the form 

tan fed =— - — —± — . 11 

E(E -V)- uuk 2 x v ; 

The dependence of energy on A vanishes, therefore, there is no parity splitting. Equation (11) 
reduces to the most simple form at V = and u = u 

, , VA 2 - u 2 k 2 . . 

t&nkd=- . 12 

uk 

From the last equation, one can simply find the number n e (h) of branches in the electron (hole) 

size-quantization spectrum (in given case n e = nu = n). This number is to be to valid an 
inequality 5 

n(n - 1) < — < irn. (13) 

It is evident from the graphical solution of equation (12) that the theorem of nonrelativis- 
tic quantum mechanics remains valid which states that there exists at least one level in an 
arbitrarily shallow ID symmetric quantum well 6 . 



5 It is interesting to note that there exists the solution to equation (12) in the case ^ = irn. This corresponds 
to a level lying on boundary with the continuous spectrum. 

6 In the nonsymmetric quantum well as in a nonrelativistic case (see, e.g., [24]), there can exist no one level. 
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The nonsymmetric quantum well of the first type. 

So we refer to the quantum well with A x = A 3 = A, A 2 = 0, V\ = —V 3 = V, V 2 = 0, u\ = 
u 3 = u, u 2 = u. Then, A explicitly enters into equation (10) as Xk x . Hence, parity splitting 
(spin splitting) arises. Energy extrema are not at k x = and lie at k x = Xk*. Hence, an energy 
depends on k x — Xk*. We suppose that \V\ <C A. In the lowest order of V, the dependence of 
energy on a momentum k x is given by expression (8) 



= ±u^m 2 + k 2 , (14) 
where k^ is the root of equation (11) at V — 



(0) _ gyg^ggg+g^gg 

tanA; d ~ ukm + (u-u)e x • (15) 



In the next order of V, we have 7 



= ±usjm 2 + (k x - Xk*) 2 . (16) 

With the same accuracy, the extremum energy ±m/c < -°- ) at k x = Xk* coincides with the extremum 
energy at k x = in the symmetric quantum well, and a difference between ones 6E is a 
higher-order infinitesimal, namely, \5E\ u\k*\ A. Expanding (15) to the second-order 
infinitesimal and remaining terms with k* 2 and k*V as terms of the same infinitesimal, we 
obtain 



u b 




(17) 
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where F = 2^A 2 - u 2 k 2 + -^=. 

The nonsymmetric quantum well of the second type. 

So we refer to the quantum well with A x ^ A 3 , V\ — V 3 — V, V 2 — 0, A 2 = 0, Ui = 
% = u, u 2 = u. We will suppose that a deviation from the symmetric quantum well is small 
Ai — A3 = 2A', |A'| A13. It is convenient to introduce the half-sum of the energy gaps 
A = Al ^ A3 , and the energy gaps are expressed in terms of A and A', Ai = A + A', A3 = A — A'. 

Analogously to the first type of the quantum well, we find the approximate solution in the 
form 

E^ = ±u\]m 2 + (jfc~ - A A;*) 2 , (18) 

where k^> is the root of equation for the symmetric quantum well (11), in which A is taken 
instead of A. Further, by analogy with previous case, one can obtain the formula 

u a 

7 The quadratic dependence of an energy on a momentum, which is valid at large values of A 2 , was used in 
Ref. [22]. 
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where a = 2(A + uk - V) \ A 2 - (uk - V) 2 + uk (uk - V)J &=((«- u)k + 

V)(A + w/c — V), c = ukniukn — V)— — , 2A "^"^l =; fcn is the root of equation (11) at 

7V " V ' (A-uk +V)y/A2-(uk -V)2 

k x = (we neglect the dependence of k^ on A^). 
T/ie interface states. 

The interface states (IS) were first predicted by I. E. Tamm in 1932 [25]. This states are 
localized near the narrow-gap semiconductor-graphene interfaces 8 . The solution in the ranges 
z < and z > d is found in the same form as for the size-quantization states and in the 
quantum well in the form 

il> xl = B ie - q * + ij X3 = B 3 e-« z + S 3 e«*. 

The quantity q is connected with E\ by the relation 

El = u\kl - u\q\ (20) 

therefore, the condition of an IS existence is the inequality 9 

\q\ < \k x \. (21) 

The equation in q is analogous to (10) with only difference that there is tanhqd instead of 
tan kd and q enters into A instead of k. 

In order to valid inequality (21) for all values of k x , it is necessary that the solution q to 
this equation vanishes at k x — > 0. Otherwise, the IS arise at some critical k x or they are not 
exist. It is simply shown that there is no such point k xQ ^ that ^(fcxo)] = \k x a\ in the case 
q{k x ) 7^ const. The function graph |?(/cx)| at k x ^ lyes either lower the straight line q = \k x \ 
(there exist the IS) or higher the line q = \k x \ (there exist no the IS). The analysis of the 
solution of the equation at small k x is presented in Table 1. 

In particular, the existence of the IS of one parity is possible. A consequence of this fact is 
a spin current along the interface. The mentioned effect can be used in spintronics [27]. 

The IS can exist for charge carriers of one sign. The IS are energetically more favorable 
than the size-quantization states. This fact results in accumulation of the IS by the charge 
carriers of one type. Thus at an optical excitation, a charge of one sign will be localized along 
the interface against the background of the uniformly distributed charge of opposite sign. 

Recently, it was shown within the electrostatic approach that the charge is induced at 
boundaries of the single GNR on the substrate at applying of an electric field transversely to 
its plane [28]. We think that a similar effect is also possible in the considered heterostructure 
due to a predominant occupancy of the IS. 



8 We note that analogous states arise in a single heteroj unction with an intersection of the dispersion curves 
[23, 26]. They are not a specificity of the quantum well and can also arise in the quantum barriers [7]. 

9 Such inequality were considered in Ref. [11] where term "evanescent modes" is used instead of the term 
"IS". 
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Table 1. The results of the analysis of the solutions of the equation for the IS with \q\ < \k. 





There are no the IS 


There are the IS 
of only one parity 


There are the IS 
of both parities 


Electrons 


j 
1 


' Pi > 1 

> 1 


1 


A = +l 
Pi>l 

,P3<1 

Pi T^PS 


A = -1 

f 

Pl<l 
[P3>1 


1 


'pi < 1 

,P3 < 1 


Holes 


1 


' Pi < 1 
k P3 < 1 


1 


A = +l 
'Pl<l 
,P3>1 

Pi 7^P3 


A = -1 
'pi>l 

,P3<1 


I 


'pi > 1 
,P3 > 1 



Notations pi = y ^-vj ' P 3 = y A3 -If are introduced. The particular case p\ — ps — 1 
corresponds to a lack of the IS at small k x for both electrons and holes (see the text). 

The case Pi = P3 = 1 (Ai = A 3 and V\ = V 3 — 0) corresponds to the lack of the IS at small 
k x for both electrons and holes. The IS of both parities can exist at \k x \ > k xc . In the case 
U\ = U2 = 113 = u, it is simply obtained from (12) 

, , \/A 2 + u 2 q 2 , x 

tanh qd = - ^ — — , 22 

uq 

the right-hand part of (22) is negative and equation (22) has not solutions. Therefore in the 
case pi = P3 = 1, an existence of the IS is possible only at it ^ « and \k x \ > k xc . 

Now, let us analyse a formation of the IS at a change of the width d of the quantum well. 
They can disappear at some critical width d c so that there are no the IS at d < d c . 

Let us determine a criterion of the existence of the IS at the defined width d. The equation 
in q has the form 

tanhx = /(x), (23) 

where x = qd and f(x) tends to an asymptotic f(x)~x at large x. For the existence of the 
solution to equation (23), the following two-sided inequality must be valid 

< /'(0) < 1. (24) 

In the general case of the nonsymmetric quantum well, inequality (24) can not be solved 
analytically and we present the results for the symmetric quantum well. 

Let us consider the case of the symmetric quantum well with V = 0. This corresponds to 
Pi — P3 — 1- The IS arise only at u > u, otherwise /'(0) < 0. One can see that the qualitative 
criterion of the IS existence, the intersection of the dispersion curves for the conduction and 
valence bands of the narrow-gap semiconductor with the dispersion straight lines for graphene 
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is realized. It should be expected that the IS exist for both electrons and holes. As it was 
noted above, the IS arise at momenta of the free motion higher some critical momentum 

1/2 



,, , 1 u + u \ L 4d 2 A 2 . 

fel><; " = 7sfc(v 1 + (^iF-^ ' (25) 

The IS exist if the inequality |g(fcar C )| < \k x \ is valid. It provides along with the condition of the 
lack of the intersection ^(A^)! with the straight line q = \k x \ a fulfilment of inequality (21) in 
the whole range I 

The fact, that the IS begin at some critical momentum of the free motion, is in accord with 
the result obtained for A x = A 2 = A 3 and V\ = V 2 = V 3 in Ref. [29]. 

For 7^0 and u — u, the IS arise at E\V < 0. They exist for electrons when V < 0, i.e. 
Pi — Ps < 1; an d for holes when V > 0, i.e. Pi = P3 > 1. There are the intersections of the 
straight lines of graphene with the dispersion curves either for the conduction band (there is 
the IS for electrons) or for the valence band of the narrow-gap semiconductor (there is the IS 
for holes). The critical momentum is 

k xc = (-u + + ( A 2 _ V 2 )d 2 ) . (26 ) 

The optical transitions. 

In presence an electromagnetic wave A = eAe 1 ^"^ where e is the polarization vector of 
the electromagnetic wave, the Dirac equation takes the form 

{ua°<y (p - ^A) + 7°Ai + Vi} * = i JU. (27) 

The operator O = — %i7°7A is considered as a perturbation. The matrix element of this 
operator is expressed by the matrix element of the rate v = e • v for the interband transition 
from the valence to conduction band E\ — ► E' x , [18] 

v = u i (* A /|7°7|* A ), (28) 

where Ui implies that U\, u 2 and m 3 are used for z < 0, for < z < d, and u 3 for z > d, 
respectively, at calculation (28). 

The calculation of the matrix element of the rate by formula (28) for the size-quantization 
gives 



p 2ikd i i 2ikd 

+f*btb)^— (fb* + b + fb*bS 



2ik w + ' 2ik 



+ i(c*c-cc*)^e- 2k * d \, (29) 



v z = 25 XX 'A* 1 A 1 | (a + a*)^- + u 2 [(b*b + + bb\ - f*f(b*_b + b^b*))d + (f*b*_b- 
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p 2ikd i i p —2ikd~ 

-rVb + )——— + (/6*6_ - fblb) 



2ik y + ' 2ik 

The coefficient A\ is determined by normalizing condition (9) 

1 + a* a „ „ „, ^ ~~ , 6*6 



+ (c*c + ?c)^-e- 2fe3 4. (30) 



2A\A X = [ + 6*6(1 + + + ^f(r( fo2 - 1) + /(( & *) 2 - !)) sin(2A;rf) + 

(r((&*) 2 - 1) - f((b*) 2 - 1)) sin 2 (fcd) + ^pe- 2 ^" , 



where 



-wifei + iXuxK _ ~ _ ~ _ fu[ 1 ~ _ u 2 k + i\u 2 k x 
E x + A 1 -V 1 '+-'-- ' ~V« 2 1-/' " 



-e , c 



/ = 1 + 



« 3 1 - / E x + A 3 -V 3 

2u 2 k(E x + Ax - K) 



—u 2 k + i\u 2 k x )(Ex + Ai — V\) + E\(iuik\ — i\uik x ) 
For the calculation of the IS matrix element of the rate by formula (28), one should replace 
k with iq 



c.,. = 2A5 AA /^i^i <{ i(a* - a)^- + iw 2 



-66+ + 2/6(6+ - b)d + fbb. 

q 1 q 



+iu 3 c(c* - c) ^ j , (31) 

^ = 0. (32) 
Here, normalizing (9) results in the expression 

2AJAi= ^-^ + 2/(6 + 6_-6 2 )d+(6 2 + 6 + 6+)^ + / 2 (6 2 + 6*6_) 



2ki v + ' + 7 2q J y ~ ' 2q 

p -2k 3 d\ -i 

+ (2 2 + C* C ) 



2fc 3 

It is evident from formulae (29)- (31) that only the transitions with the conservation of the 
parity is allowed (Fig. 2 shows them by arrows). The optical transitions for the IS exist only in 
the particular case p\ = p 3 = 1 at an absorption of the electromagnetic waves with the linear 
polarization along the interfaces. 

These conclusions can be experimentally verified by the following way. Let us take the 
symmetric quantum well. Applying a voltage U to the GNR, one can achieve that the work 
function V becomes equal to zero, i.e. the particular case pi = p 3 = 1 is realized at U ^ 0. 
Then the optical transitions arise. They are due to the absorption of the electromagnetic 
waves with an energy E so that E min < E < A e y/ where A e ff is the effective energy gap in 
the size-quantization spectrum (the smallest energy difference between the electron and hole 
branches). 
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The possibility of the exciton formation. 

In conclusion, we note that a formation of the finite effective energy gap in the size- 
quantization spectrum can result in the exciton generation at the optical pumping (the energy 
gap in graphene is zero and there are no excitons). 

The ortoexcitons with the total spin S = 1 of an electron and a hole are direct in the the 
symmetric quantum well therefore they are short-lived, and the paraexcitons with S = [30] 
are long-lived, since the transitions with the change in the parity are forbidden (out of the 
dependence wether they are direct or not). The ortoexcitons are indirect in the nonsymmetric 
quantum well (Fig. 2), therefore they must also be long-lived. It was noted earlier that the 
paraexcitons are short-lived for the layered narrow-gap semiconductor heterostructures because 
the transitions with the change of the parity are allowed [31]. 

The calculation of exciton binding energy in the planar quantum well based on graphene and 
the narrow-gap semiconductors represents a separate problem and it will be present elsewhere. 




Fig. 2. Splitting of the energy spectrum in the nonsymmetric quantum well. The allowed 
transitions with the conservation of the parity are shown. 
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